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We investigate the entanglement patterns of photon-added and -subtracted four-mode squeezed 
vacuum states. Entanglements in different scenarios are analyzed by varying the number of photons 
added or subtracted in certain modes, which are referred to as the “player” modes, the others 
being “spectators”. We find that the photon-subtracted state can give us higher entanglement 
than the photon-added state which is in contrast of the two-mode situation. We also study the 
logarithmic negativity of the two-mode reduced density matrix obtained from the four-mode state 
which again shows that the state after photon subtraction can possess higher entanglement than that 
of the photon-added state, and we then compare it to that of the two-mode squeezed vacuum state. 
Moreover, we examine the non-Gaussianity of the photon-added and -subtracted states to find that 
the rich features provided by entanglement cannot be captured by the measure of non-classicality. 


I. INTRODUCTION 

Distribution of entanglement in a multipartite quan¬ 
tum system is known to be a useful resource in several 
quantum communication and quantum computational 
tasks [1]. Notable ones include quantum secret sharing 
[2], distributed quantum dense coding [3], distribution 
and concentration of quantum state [4], and cluster state 
quantum computing [5] . Such protocols have successfully 
been realized in physical systems like photons [6], ions [7], 
nuclear magnetic resonance [8] , nitrogen vacancy centers 
[9], etc. 

One of the physical systems in which quantum informa¬ 
tion tasks have been realized in the laboratory is the class 
of continuous variable (CV) systems. Historically, the no¬ 
tion of the quantum correlated state of two particles in 
CV systems first arrived in the seminal paper of Einstein, 
Podolosky, and Rosen in 1935 [10]. In recent years, sev¬ 
eral communications schemes like teleportation [11] and 
classical information transfer by quantum channels [12], 
have extensively been investigated both theoretically and 
experimentally, in CV systems, especially in Gaussian 
states [13-16]. However, it has been discovered that there 
are several protocols which can not be implemented using 
Gaussian states with Gaussian operations. Examples in¬ 
clude entanglement distillation [17], measurement-based 
universal quantum computation [18], teleportation [19], 
and quantum error correction [20] . The increasing impor¬ 
tance of non-Gaussian states have led to the discovery of 
several mechanisms to create such states in the labora¬ 
tory [21], An important one is adding and subtracting 
photons, when the initial state is the squeezed vacuum 
state. Starting with the single mode squeezed vacuum 
state, whose Wigner function [22] is always positive, it 
was shown that photon addition can generate a negative 
dip of the Wigner function in the phase space [23] and 
hence can deviate from being a Gaussian state. In case 
of the two-mode squeezed vacuum (TMSV) state as in¬ 
put state, both entanglement and fidelity of teleportation 
can be increased by adding and subtracting photons to 
(from) one or two modes [24]. For such experiments, see 


[25]. Moreover, the entanglement content of the photon- 
added state obtained from the TMSV state was shown 
to be always higher than that of the photon-subtracted 
state [24]. 

Investigations of the squeezed vacuum state with re¬ 
spect to photon addition and subtraction are usually re¬ 
stricted to the two-mode case, even though the impor¬ 
tance of multimode CV system is unquestionable. In 
this paper, we consider the four-mode squeezed vacuum 
(FMSV) state as input and deGaussify it by adding and 
subtracting photons in different modes. We evaluate en¬ 
tanglement between different modes in all possible bi¬ 
partitions and compare the results of the photon-added 
state with the subtracted ones. We call a mode as “spec¬ 
tator” mode in which either no photon or fixed number 
of photons are added (subtracted). The other modes are 
referred to as “players”. We here investigate two sce¬ 
narios - (1) one player mode (2) two player modes. We 
analytically show that in the single player case i.e., when 
photons have been added (subtracted) to (from) a sin¬ 
gle mode, in the player : spectator bipartition, entangle¬ 
ments of the photon-added and -subtracted states coin¬ 
cide. In this situation, we prove that entanglements in 
both photon-added and -subtracted states monotonically 
increase with the number of photons added or subtracted. 
Unlike the TMSV case, we observe that there exists sce¬ 
narios in which photon-subtracted output states obtained 
by subtracting photons from one or two modes con¬ 
tain higher entanglement, compared to the photon-added 
state. Specifically, we find that the photon-subtracted 
state contains more entanglement in the spectator : rest 
bipartition than that of the photon-added state, when 
a single mode acts as a player. Similar hierarchy can 
also be obtained when any two modes act as players. 
Interestingly, the advantageous situation for photon sub¬ 
traction can be reversed by adding fixed number of pho¬ 
tons in the spectator modes. Such behavior can also be 
viewed by analyzing logarithmic negativity of the output 
two party state which can be obtained by discarding ei¬ 
ther two player modes or two spectator modes. Finally, 
we study a distance-based measure of non-Gaussianity 
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in these scenarios and find that the non-Gaussianity in 
general is higher for the photon-added state than that 
of the photon-subtracted state. In case of two inodes, 
photon-added states are known to be more non-Gaussian 
than the photon-subtracted states. However, as already 
noticed in [24], the relation between entanglement and 
non-Gaussianity is not straightforward. In the four-mode 
case, we again find that the photon-added state has al¬ 
ways higher non-Gaussianity than that of the photon- 
subtracted states, and hence reflects that entanglement 
and non-Gaussianity are possibly not directly connected. 

The paper is organized in the following way. In Sec. 
II, we discuss briefly the N-mode squeezed vacuum state, 
and two special cases, the two mode squeezed vacuum 
state, in II A, and the four-mode squeezed vacuum state, 
in IIB. In Sec. IIB 1, we consider the FMSV state, 
when TOj number of photons are added to or subtracted 
from the mode i. In Secs. Ill and IV, we briefly in¬ 
troduce a non-classicality measure of a quantum state 
in a continuous variable (CV) system, and the quan¬ 
tum correlation measures which are relevant in the pa¬ 
per, respectively. In Sec. V, we present the main re¬ 
sults in which we systematically compare the entangle¬ 
ment of the four-mode photon-added state with that 
of the photon-subtracted state by considering the von 
Neumann entropy in different bipartitions. Another en¬ 
tanglement measure, the logarithmic negativity, for the 
photon-added and -subtracted states are evaluated in 
Sec. VI, while the behavior of non-Gaussianity for the 
output state is studied in Sec. VII. We summarize in Sec. 
VIII. 


II. N-MODE SQUEEZED VACUUM STATE 


In this section, we discuss the N-mode squeezed vac¬ 
uum state (NMSV), specifically the two-mode and four¬ 
mode squeezed vacuum states, and a state obtained af¬ 
ter adding (subtracting) an arbitrary number of photons 
at the mode i. These states are examples of entangled 
states in continuous variables which can be used in var¬ 
ious quantum information tasks. To define such states, 
let us first denote the bosonic creation and annihilation 
operators at the mode i, as dj and di respectively, which 
satisfy the bosonic commutation relations, [dj,dj] = 5-ij. 

and [d,;, dj] =0, [dj, dj] = 0. By using bosonic operators, 
an TV-mode squeezing operator can be defined as 


v i— 1 


e*difii + i - edjdj +1 


where djv+i = di. The corresponding (NMSV) state is 
given by 

IV’jv) =«S(e)|0 1 0 2 ...0 JV ) 

1 ( i N 

= T\T eXP 1 ~2 ^ dj tanh(rQ) ife djy e 

\ i. fe =i 

|0i0 2 ...(hv), (2) 

where |0i0 2 ... On) is the N-mode vacuum state, Ns, is 
a normalization constant, and e = re , with r being the 
squeezing parameter. Here the matrix Q is obtained from 
the following relation 

N 

«S(e) t d J <S(e) = ^(cosh(rQ) jfc d fe - sinh (rQ) jk e l9 a[) 

k 

V j = l,...,TV(3) 

Let us now define the position and momentum operators 
for each mode, given by 

Qi = (d, + dj), (4) 

Pi = v(dj — dj), (5) 

% 

to show that the Eq. (2) indeed represents a squeezed 
state. The variances of the N-mode quadrature opera¬ 
tor 8 ’ and *2 = 2^57 Ej^-aj), 

are given by 

AXf = i(e 2r sin 2 (0/2) + e -2r cos 2 (0/2)), (6) 

and 

A AT 2 = ~(e 2r cos 2 (0/2) + e _2r sin 2 (0/2)). (7) 

Thus for 9 = 0 or it, we have AATiAAT 2 = ). However, 
for anyone of the i = 1,2, A Ad < j, while A Ad > \ for 
the other i. This guarantees that the state, given in Eq. 
(2) is a squeezed state. We assume 9 = 0 throughout the 
paper. 


A. Two mode squeezed vacuum state 

The two mode squeezed vacuum state can be obtained 
by putting TV = 2 in Eq. (2), where Q = ^ d J ^ , and 

tanh ^ = ( tanh r ^ ‘ 

Thus, the TMSV state with 9 = 0 is given by 
|V> 2 > =sechre- tanh ™^|00) 

OO 

= sech r ^~j(— tanhr) n |re)|n), 

n=0 


S(e) = exp 


( 1 ) 


(9) 
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where | n) = | 0 ), is the occupation number state. 

Taking \ip 2 ) as the initial state, the behavior of entan¬ 
glement and non-Gaussianity after adding or subtracting 
photons, have extensively been investigated [24], 


B. Four-mode squeezed vacuum state 


Expanding the exponential in Eq. (11), we have 


IV’t) = 


cosh; 


£ - 


n —0 



Let us now consider the FMSV state obtained by set¬ 
ting N = 4 in Eq. (2). The 4x4 matrix, Q , in this case, 
takes the form 


Q 


1 

2 


/ 0 1 0 1 \ 
1 0 1 0 I 
0101 
V1 0 1 0 / 


( 10 ) 


The FMSV state with 9 = 0, is then given by [26] 


l^4> = 


1 


coshr 


: (al+4)(4+4)|0000). ( 11 ) 


1. Photon-added, and -subtracted four-mode state 


In this paper, we consider the FMSV state, ^ 4 ), as 
an initial state and our aim is to find the characteristics 
of its entanglement and the measure of non-Gaussianity 
after adding and subtracting a finite number of photons. 
Suppose rrii number of photons are added at each mode 
i, with i = 1,2, 3,4. Then the output four-mode (FM) 
state reads as 


n= 0 ' ' r\ ,r2— 0 V 


n\ /n\ /(n — ri + mi)! /{n - r 2 + m 2 )! /(ri+m 3 )! /(r 2 + m 4 )! 


(n - n)\ 


(n - r 2 )! 


r 1 ! 


rp- 


00 n 


= E Pn™i!r 2 l n - n + mi)|n - r 2 + m 2 )|ri + m 3 )|r 2 + m 4 ) 

n —0 ri ,r2—0 


x | n — n + mi)|n — r 2 + m 2 )|ri + m 3 )|r 2 + m 4 ) 

(13) 


where N add is the normalization constant. Similarly, after subtracting {rrii}(i = 1,2, 3,4) number of photons from 
each mode of the FMSV state, the resulting state is given by 


cxj / . , \ n n—m i n—m 2 

1 «>> = £ (-^) E E 


n—M 


00 n—m 1 n—m 2 


n —m 3 r 2 —m .4 


(n - r 4 )! 


(n — r 2 )! 


n! 


r 2 ! 


fiy \ r 2/ V (« — r 4 — m 4 )! V (n — r 2 — m 2 )! V (r 4 — m 3 )! V (r 2 — m 4 )! 

x |n - r 4 - mi) |n - r 2 - m 2 ) |ri - m 3 ) |r 2 - m 4 ) 


= E E 1 71 -’"l-mi)|n-r 2 -m 2 )|ri-m 3 )|r 2 -m 4 ), 

n—M r 1=777,3 r-2—m 4 


(14) 


where 7V sub is the normalization constant, and M = 
maxjmi + m 3 , m 2 + m 4 }. 

III. MEASURE OF NON-CLASSICALITY IN 
CONTINUOUS VARIABLE SYSTEMS 

The negative Wigner function of a given state indicates 
the non-classical nature of the corresponding state while 
the positivity implies the opposite. On the other hand, it 
is known that the Wigner function of a Gaussian state is 
always positive [27]. Therefore, one can define a measure 
of non-Gaussianity or non-classicality by measuring the 
departure of a given state, p , in a CV system from a 


Gaussian state. In terms of relative entropy distance, it 
is given by [28-30] 

$ ng ( q ) = ^(elleG) 

= S(qg)-S(q), (15) 

where Seller) = —tr( 7 ylog 2 a) — S(r 7 ), and pa is a Gaus¬ 
sian state which has same covariance matrix and first 
moment as p. Here, S(a) = — tr(crlog 2 ( 7 ) is the von 
Neumann entropy of er. 

The von Neumann entropy, S(qg), of any Gaussian 
state can be calculated by using its covariance matrix, a. 
For an N mode Gaussian state, qg , the von Neumann 
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entropy is defined [14] as 

N 

S(ea) = ^2g(^k), (16) 

k= i 

where Vk is the Williamson normal form of the covari¬ 
ance matrix of the iV-mode Gaussian state go, and the 
function g(x) is given by 

d ( x ) = 2 lo ® 2 ( 2 / 2 °^ 2 ( 2 / (17) 

In this paper, we will take a Gaussian state as the 
input state, and after photon addition (subtraction), the 
diversion of the output state from the input Gaussian 
state will be quantified by Sng- 


IV. QUANTUM CORRELATION MEASURES 

Quantum correlation measures in bipartite systems, es¬ 
pecially for two qubit systems, are well understood. Such 
quantifications include the von Neumann entropy of lo¬ 
cal density matrices for pure states [31], entanglement 
of formation [32], concurrence [33], logarithmic negativ¬ 
ity [34]. However, measures of quantum correlation in a 
multipartite scenario, both in discrete and CV systems 
are limited [1, 15]. To characterize entanglement in CV 
system with multiple modes, one possibility is to compute 
von-Neumann entropy in different bipartition of modes. 
Another possibility is to study logarithmic negativity of 
two modes which can be obtained after discarding all the 
modes except two. In this section, we briefly discuss the 
local von Neumann entropy and the logarithmic negativ¬ 
ity in CV systems. 


is in a pure state. The entanglement between two blocks 
of a multiparty state can capture entanglement distribu¬ 
tion in the multipartite domain. Such a quantification 
has been extensively used in many-body systems [35]. 
Here, we divide the multimode system into two parts 
and investigate the behavior of the entanglement con¬ 
tent in the bipartition by adding (subtracting) photons 
in various modes. 


B. Logarithmic negativity 


In CV systems, logarithmic negativity (LN) is an im¬ 
portant entanglement measure [34]. For a state, pn, with 
N = Ni + N 2 modes, it is given by 

LN( Pn ) = log 2 J\f(p N ), (19) 

where the negativity of the given state is given by 


N(pn ) — 1 + 2 Pi 


( 20 ) 


Here /q’s are the negative eigenvalues of the partially 
transposed density matrix, p ^ 1 , where partial transpo¬ 
sition is taken with respect to the N\ modes [36]. As 
mentioned for the evaluation of the von Neumann en¬ 
tropy, LN is also calculated by truncating to a large block 
of the infinite dimensional matrix. 


V. COMPARISON OF ENTANGLEMENT 
ENHANCEMENT BETWEEN PHOTON 
ADDITION AND SUBTRACTION 


A. Entanglement of a pure state 

Entanglement of a bipartite pure state, can be defined 
by the von-Neumann entropy of the reduced density ma¬ 
trix of a given state [31], i.e., 

E{ W)ab) = S(p A ) (18) 

where pa = tred^AB^I)- In CV systems, entangle¬ 
ment of a two-mode pure state can be quantified by 
the von Neumann entropy of a single mode. The sin¬ 
gle mode density matrix can be a matrix of infinite di¬ 
mension which has to be diagonalized to evaluate its von 
Neumann entropy. The calculation of the entropy can be 
carried out after truncating the matrix to a large block. 
The block size is determined by checking for convergence 
of trace, with increasing block-size, to unity up to a cer¬ 
tain significant digit. We will discuss this issue in detail 
for a specific scenario. 

In the multipartite domain, entanglement is difficult to 
characterize even for pure states [1]. However, if one di¬ 
vides a multipartite system into two blocks, then the en¬ 
tanglement between the two subsystems is the von Neu¬ 
mann entropy of one of the block, provided the system 


In this section, our aim is to investigate the effects on 
entanglement in different bipartitions, when photons are 
added (subtracted) in (from) different modes of a four¬ 
mode squeezed vacuum state. To study such behavior, 
we divide the modes into two different categories, viz. (1) 
player modes - the modes in which number of photons 
that we add (subtract) varies, and (2) spectator modes 
- the modes in which either no photon or fixed num¬ 
ber of photons are added (subtracted) and hence plays a 
spectator role in the deGaussification process. The com¬ 
parison has been made between the situations, when the 
rrii, i = 1,2 ,... photons are added in the player modes, 
and the scenario when the same number of photons are 
subtracted from the player modes. To execute such com¬ 
parison, we introduce a quantity 

«5(M) = EffC 1 - 1 ) - BKiJ" 1 ) (21) 

where A: B is a bipartition with A D B = 0. The pos¬ 
itivity of S E ({rrii}) implies that addition is better than 
subtraction from an entanglement perspective. It is clear 
that the behavior of ^({uii}) with {to,} depends on the 
number of player and spectator modes as well as the bi¬ 
partite splits. 
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FIG. 1. (Color Online) Schematic diagram of choices of player 
and spectator modes as well as partitions. If we fix the bi¬ 
partition to be 1 : 234 there are three nontrivial possibilities 
of choosing a single player in the photon-added and the sub¬ 
tracted FM state. There are the cases (a) - (c), and the 
number in the square mentioned for each case is the mode at 
which the photon is added/subtracted. 


A. Photon added and subtracted with one player 
mode 


Let us first consider a situation in which one mode acts 
as a player while the rest are the spectator modes. We 
first restrict ourselves in the 1 : 234 cut irrespective of 
the choice of the player mode. In this case, there exists 
three different possibilities of choosing a player mode (a) 
first mode as player and the rest as spectators, ( b ) second 
mode as player and (c) third as player (see Fig. 1). From 
Eqs. (13) and (14), it is clear that fourth mode as a 
player is equivalent with case ( 6 ), and hence we exclude 
this case. 


1. Single player mode in the smallest bipartition 


Suppose that we add or subtract mi photons in the 
first mode without putting any number of photons in 
rest of the modes, as shown in Fig. 1(a). Here the first 
mode acts as a player. The reduced density matrices can 
be calculated from Eqs. (13) and (14), which read as 


add 1 O tanh 2T V O /n\ (re + mi -Tq)! 

i' mi 2 " (n-ri)! 

\n + mi — T\)(n mi — r\ | (22) 


for photon addition, and 

sub _ 1 y-tanh 2n r ^ (mi+ri)!/n + mi\ 

Pi.rm - N s ub 2 " ^ rq! Ui+mW 

1 n=0 ri—0 N ' 

l r i)( r i| (23) 

for photon subtraction. We now analytically establish 
that entanglement in the bipartition of the player and the 
spectator modes, increases with the number of photons 
added. 

Proposition: Entanglement increases with the ad¬ 

dition of a single photon in a four-mode photon-added 
state , i.e., 

£(I< + 1»1:234 > E(|< d »i :234 (24) 

where = 0, 1 denotes the state in which m\ + i 

number of photons are added at the mode 1 . 

Proof : To evaluate entanglement in the 1 : 234 bipar¬ 
tition, we have to study the single mode reduced density 
matrix, pi rl r d ni , of the four-mode state |Off )■ To prove 
E(\ip^ d +1 )) 1:234 > £( | of)) 1 : 234 , it is equivalent to show 
^(Pifmi+i) — After inserting the normaliza¬ 

tion constant in Eq. (22), we get 


Pi, 


add _ 2 mi 


i 1 -*)* 1 * 1 f „ , 

(2-sp. ,0 ( 


m i + r i 

mi 


|mi + ri)(mi + rq| 

OO 

= ^ 3 (x,mi,ri)|mi+ri)(mi+ri| (25) 

rq —0 


where x = tanh 2 r, 


f(r,x) 



and 


g(x, m,r) = 2 m 


(1 -x ) m+1 
(2 — x) m 


f{r,x) 


fm + r\ 

V m )‘ 


(26) 


(27) 


Therefore, entanglement in the player : spectator bi¬ 
partition is given by £(|Of)i:re*t) = = 

- E^=o 9(x, mi, ri) log 2 g(x, mi,ri). 

Now if we add one more photon to the state in Eq. 
(22), the entanglement is going to be E(\ip%ff +1 )i :rest ) = 

- E"=o 9 ( x , mi + 1, ri) log 2 g(x, mi +1, n). Let us now 
evaluate g{x,mi + l,ri). It simplifies as 


g{x,mi + l,n) 


ff(x,mi,n) 


2(1-x) 

( 2 ~x) 

+ 77 — g( x ,mi + l,ri - 1 ), (28) 
2 — x 


by using Pascal’s identity, and the recursion relation of 
f(r,x), which is given by 


f(r,x) 


X 


2 — x 


f(r - 1,0 


(29) 
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Using the concavity of the function h[x) = —x log 2 x, we 
get 


2(1 — x) 

h(g(x,mi + l,ri)) > —- r-h(g(x, mi, ri)) + 

(2 — x) 

h(g(x,m 1 + l,ri - 1)).(30) 


2 — x 


Taking the sum over r\ in both sides, we have 


S(p 


add 

1 , 771 - 1+1 


)> 


2(1 -x) 
(2 — x) 


S{p 


add 

l,mi 


)+ 


X 

2 — x 


S(P^ 1+ 1) (31) 


which immediately implies 


Of „add \ \ o/ add \ 

d[Pl,m 1 + l) > M/ 3 l,mj- 


(32) 


Hence the proof. ■ 

Similarly one can also show that entanglement of the 
photon-subtracted state in the player : spectator split 
increases with number of photons subtracted from the 
state. 

We are now going to analyze the effects on entangle¬ 
ment under addition and subtraction of same number of 
photons . 

Proposition: When a single mode acts as a player, 
entanglement between the player and the spectator modes 
of the photon-added state coincide with that of the 
photon-subtracted state. 

Proof : To prove that the increase of entanglement in 
the multimode state is same for addition and subtraction, 
we consider the single mode reduced density matrix. The 
single site reduced density matrix of photon-subtracted 
state after inserting N-f ub , is given by 


Pi, 


sub _ _ jjmi+1 \ ' 


7 * 1=0 

oo 

<£ 


mi + ri 


n + mi 

ri + mi 


|ri)(n|,(33) 





FIG. 2. (Color online) Behavior of E( \ipm d ) 1.234) and 
E( IVC+ , 234 ) vs. mi. We add (x) and subtract (+) upto 
40 photons in (from) the first mode, and calculate entangle¬ 
ment in the 1 : 234 bipartition, when no photons are added 
(subtracted) in (from) the spectator modes. As shown in the 
propositions, entanglement in both the cases increases mono- 
tonically with mi and they coincide. 


To visualize the above Propositions, we plot 
S(pi d ^/ 1 sub ), with respect to mi by fixing the squeezing 
parameter r = 0.4 in Fig. 2. It clearly shows that the 
curve for photon addition merges with the curve of pho¬ 
ton subtraction. Moreover, it shows that entanglement 
in that bipartition monotonically increases with the ad¬ 
dition or subtraction of photons as shown in Proposition 
1. Note here that although the results presented here 
are when the photons are added at the mode 1 and the 
bipartition is considered as player : spectator mode, the 
Propositions remain unaltered if another mode also acts 
as a player by keeping the similar bipartition. 


2. Ejfects on entanglement due to change of partition 


where 

2 m i+i 

fsub (ri, mi, x) = (2 1 x y 1+ m 1+ i > ( 34 ) 

which can be obtained by a recursion relation similar 
to that given in Eq. (29). On the other hand, the re¬ 
duced density matrix after adding same number of pho¬ 
tons reads as 

add y nmi (l^) mi+1 2 / * yi 

Pl ' mi ^ (2- 2+1 2 -x\2-x) 

X ( Wl m~ ri ) l ri + mi )( ri + TOl l- ( 35 ) 

Comparing Eqs. (33) and (35), we have S(pi d f} li ) = 

S{p{%). ■ 


We now consider the entanglement in the same bipar¬ 
tition as in the previous case, i.e., 1 : 234. However, 
the second or third mode now act as player and no pho¬ 
tons are added in the rest of the modes. In the previous 
case, one block contained only the player mode while the 
other one contains all the spectator modes. In this case, 
one part of the partition contains one spectator mode 
while the other one consists of both the player and the 
rest of the spectator modes. In the previous case, we 
have already shown that the effects on entanglement due 
to addition and subtraction of photons are similar. We 
will now show whether such observation remains invari¬ 
ant even in this scenario. 

Let us now take the four-mode squeezed vacuum state 
as input, and add (subtract) m 2 photons in (from) the 
second mode. As depicted in Fig. 3(a), we find that un¬ 
like the previous case, the photon-subtracted state pos- 
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m 2 m 3 


FIG. 3. (Color online) (a) Trends of E{\ip^ d ) 1 : 234 ) and 
£(|VC 2 6 >1:234) with the number of photon-added (subtracted) 
in (from) the second mode, (b) Similar study has been car¬ 
ried out when the third mode acts as a player. Both the cases 
reveal that subtraction is better than addition. 


sesses more entanglement in the 1 : 234 bipartition than 
that of the photon-added state. The ordering remains 
unchanged if one takes the third mode as player and con¬ 
sider entanglement in the 1 : 234 split (see Fig. 3(b)). 
Moreover, we observe that the amount of entanglement 
decreases in this scenario, compared to the case when 
the second mode acts as a player. Note here that if one 
takes the two-mode squeezed vacuum state as input, it 
was observed that the bipartite entanglement content of 
the photon-subtracted state is always lower than that of 
the photon-added state. 

3. Bipartition with both player and spectator modes 

We still restrict ourselves to the case of a single player. 
But we now move to the situations in which entangle¬ 
ment of a four-mode state is studied by considering a 


bipartition in which both sides of the split contain two 
modes, namely 12 : 34 and 13 : 24. The other split be- 



FIG. 4. (Color online) Schematic diagram of two different 
blocks, when a single mode, specifically the first mode, acts 
as a player. 

tween modes, i.e., 14 : 23, reflects a similar behavior, due 
to the symmetry of the four-mode state. In these two 
scenarios, photons are added or subtracted in the first 
mode, as shown in Fig. 4, and no photons are added or 
subtracted, in the other spectator modes. 

To study entanglement of | tp^f) (l^m, 6 )) in the 12 : 34 
or 13 : 24 bipartition, we require the two party reduced 
density matrices, p ^ d mi , , p s ^ b mi , and p\f mi . We 

have 


^ ""W /FTT" 

2 ” + "' „„„ V wU 


and 


1 OO (n+nl/2 min{n,n'}- m imin{n,n'} 

sub _ 1 \ ' X \ " \ ' 

rl2,mi jysub / v 2 n + n/ / -j / -j 

n,n'—mi ri—0 V 2—0 


Note that in the previous cases, where one partition 
contains only a single mode, we required single-site den¬ 
sity matrices to calculate the entanglement, and they 
are always diagonal in the number basis. The same is 
not the case for two-site density matrices. Similarly, one 
can find out the reduced density matrices of Pi dd mi and 
Pi 3 mi- In both the scenarios, we observe that entan¬ 
glement increases against the number of photons added, 


’ (n + mi — n)! / (n' + mi — rq)! 


GT/ \ r 2y 


(n - ri)! 


(n' — ri)! 


x|n + mi - ri)i\n - r 2 ) 2 (n + toi — ri|i(n' - r 2 12 , (36) 


n 

n 



x|n — mi — ri)i|n — 


I [n — ri)\ / in’— ri)\ 

(n — mi — ri)\ V (n' — mi — n)! 

T’ 2 ) 2 (n' - mi - ri\i(n' - r 2 1 2 . (37) 


mi and same is true for subtraction of photons (see Fig. 
5). Moreover, as observed in the previous case with 
the smallest partition consisting of the spectator mode, 
photon-subtracted state contains higher entanglement in 
the 12 : 34 as well as 13 : 24 partitions than that of the 
corresponding photon-added state. See Figs. 5(a) and 
5(b). 

We briefly mention here the method used to calculate 












































FIG. 5. (Color online) Plots of entanglements of photon- 
added and -subtracted states in the 12 : 34 (a) and 13 : 24 
(b) bipartitions with mi. 


S(Pi 2 mi), an d the other local entropies. The von Neu¬ 
mann entropy of Pi dd mi can be obtained if one can di¬ 
agonalize the infinite dimensional matrix, given in Eq. 
(36). To calculate it, for fixed mi, we have to truncate 
the summation upto a large value of n and n', say N for 
both, and calculate its trace, i.e., trjv(p “2 mi)> as we U as 
von Neumann entropy, S N (pi dld mi ). We then choose, 27V 
as maximum of n and n' and obtain the quantities. When 
the difference between S N (p1 dd mi ) and S 2N {p \ dd mi ) is °f 
the order of 10~ 6 , we take S N (pi dd mi ) as the actual en¬ 
tropy. In Fig 6, for a fix value of mi, we plot S N (pf dd ) 
and tr 7 v(pi 2 <1 ) with the variation of TV. With the increase 
of mi, we observe that we require higher values of TV. 
However the figure shows both the quantities converge 
when TV > 10, irrespective of the value of mi. When we 
compute entropy or LN, we always carry out a similar 
scaling analysis for choosing TV. 



N N 


FIG. 6. (Color online) Plot of convergence of von Neumann 
entropy, S N (pi$ d mi ), in (a) and trjv^il.mj in ( b ) against N 
which is the maximum value of n and ri . We choose three 
different values of mi, viz. mi = 10, 25,40. We find that for 
example, for mi = 40, trace goes to unity and the entropy 
(entanglement) converges for N > 10. 


B. Behavior of entanglement of photon-added and 
-subtracted states with two player modes 

In this section, keeping the four-mode squeezed vac¬ 
uum state as the input state, we increase the number 
of players from one to two modes, and hence the possi¬ 


bilities of choosing the player modes with nontrivial bi¬ 
partition grows substantially. For a fixed bipartition, we 
investigate the nature of entanglement by changing the 
modes in which photons are added or subtracted. Upto 
now, we have shown that the entanglement content of the 
resulting state after subtracting photons is either equal 
or higher than that of the photon-added states. Let us 
now investigate whether such situation persist when two 
modes are players. 



FIG. 7. (Color online) Schematic diagram of four non-trivial 
possibilities of choosing two modes as players in the 1 : 234 
bipartition. Other choices can be shown as repetitions due to 
the symmetry of the FM state. 


1. One part of the bipartite split contains a single mode 

We begin by concentrating on the entanglement of the 
FM state after addition (subtraction) of photons in the 
1 : 234 bipartition. In this scenario, there are four pos¬ 
sibilities for adding and subtracting photons. As shown 
in Fig. 7, the modes that act as players are as follows: 
(a) the first and the second mode, (b) the first and the 
third mode, (c) the second and the third mode, and (d) 
the second and the fourth mode. Other possibilities can 
be reduced to any one of the above four cases due to the 
symmetry in the four-mode squeezed state. Moreover, it 
can be shown that the entanglement pattern of cases (a) 
and (b) are qualitatively similar while cases (c) and (d) 
are analogous and hence the entanglement features will 
be studied in pairs. 

Cases (a) and ( b ): We now consider the situation 
where either the first and the second modes act as players 
or the first and the third modes are players. We calculate 
the 6i(mi,mi) (i ^ 1), when no photons are added and 
subtracted from the spectator modes. We observe that 
there exists a region for which > 0, which is 

in contrast with the case when one mode was player in the 
preceding subsection (see Fig. 8(a)). As seen from the 
figure, for moderate values of mi, the boundary between 
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FIG. 8 . (Color online) Behavior of df (mi, m3) against m 1 
(horizontal axes) and m 3 (vertical axes). Panels (a) and (b) 
correspond to inactive and active spectator modes respec¬ 
tively. I 11 Fig. (b) we add (subtract) m 2 = 5 photons in 
(from) the second mode. If both second and fourth modes 
are active spectators, the region of Si(mi, m 3 ) increases in 
(mi, m, 3 )-plane. For example if we choose m 2 =, and m 4 =, 
the region for which Sf(mi, m 3 ) > 0 increases. 


(a) 



FIG. 9. (Color online) Distinct scenarios of two player modes 
in the 13 : 24 split. There are two possibilities - (a) first and 
second as players, (b) first and third as players. 


the positive and negative regions is almost a straight line 
and hence we can find the slope of the straight line which 
can help to study these cases quantitatively. We find that 
for high values of m 1 , the slope of df (m 1 , m 3 ) = 0 is ap¬ 
proximately 0.28, which is small compared to the slope 
of = 0, which is 0.64. Moreover, we notice 

that max [< 5 ^( 777 . 1 , 7773 )] = 0.2 < max [< 5 ^( 7771 , 7772 )] = 0.4 
while minimum value of fif(mi,m 3 )(= —2.0) is smaller 
than that of Sf(mi,ni 2 )(= —1-6), in the regions sur¬ 
veyed. Therefore, we can conclude that to create max¬ 
imal entanglement in this scenario, photon addition is 
advantageous when one adds photons in the first and the 
second modes compared to the case of rrii and 7773 being 
players (with 7771 rrii, i = 2,3). 

In both the cases, spectator modes play an important 
role in the behavior of entanglement in the 1 : 234 bipar¬ 
tition. As depicted in Fig. 8 (b), entanglement in the 
photon-added state can be increased by adding photons 
in the spectator modes. For example, when 7772 ( 4 ) = 5, 
<5^(7771,7773) against 7771 and 7773 is depicted in Fig. 8 (b). 
Quantitative comparison can be made between Figs. 8 
(a) and 8 (b). In particular, for 7771 3 > 7773, the region 
with < 5 ^( 7771 , 7773 ) > 0 when no photons are added (sub¬ 
tracted) in the spectator modes can be calculated. In 
this limit, we assume that the boundary is a straight line 
and hence the area is the area of a quadrilateral. Let us 
call the area as Ao- In this case, we calculate the area 
of the quadrilateral when mi > 25 and 7771 < 40, and we 
find A 0 « 160. After adding (subtracting) 5 photons in 
the second or fourth modes, we find that the area, A 5 , of 
the corresponding quadrilateral increases and A 5 « 253. 

Behavior of entanglement in the 1 : 234 split for cases 
(c) and (d) are almost identical with the previous cases. 
The only difference is that entanglement of the sub¬ 
tracted state is always better than that of the added state 
when spectator modes are inactive. The picture changes, 


i.e. entanglement of the photon-added states starts in¬ 
creasing faster than the photon-subtracted states, like in 
the preceding cases, when fixed numbers of photons are 
added (subtracted) in the spectator mode(s). 


2. Bipartition containing equal number of modes 

We will now consider the case where, we still keep two 
modes as players but we now divide four-modes into two 
blocks consisting of two modes instead of one mode in 
the preceding discussion. In this case, the two nontrivial 
bipartitions are 12 : 34 and 13 : 24. Let us first con¬ 
centrate on the bipartition 13 : 24. In this case, the 
symmetry of the FMSV state after addition or subtrac¬ 
tion of arbitrary number of photons in all the modes, 
given in Eqs. (13) and (14), ensures that there are only 
two nontrivial situations in the case of two player modes 
(see Fig. 9). They are - (a) when the players are the 
first and the second modes, and ( b ) when first and third 
modes act as players. Cases (a) and (b) show similar en¬ 
tanglement behavior like previous cases, when one part 
of the bipartition contains a single mode, and hence we 
only discuss the situation when two spectator modes are 
active, which have not been analyzed before. 

Case (a): The reduced density matrix of the first and 
the third mode, for the photon-added state, is given by 


Plt{m,}= tr 24 (|^ } )<^ } l) 

00 n 

= ^2 a n,ri,q\n + 777i - ri)i|777 3 + ri) 3 

(77 + 7771 - d|l (7773 + d| 3 , (38) 


n—0 r i=0 
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FIG. 10. (Color online) Role of spectator modes in 
<S^i(mi,m 2 ). In (a), m 3 = 10 and U 14 = 0, while in (b), 
m 3 = rriA = 5. We see that spectator modes help to enhance 
entanglement in the photon-added state. 



FIG. 11. (Color online) Schematic diagram of choices of two 
player modes in the 12 : 34 split. 


where we write 

(COG 

/(n + mi - ri)! (n + mi — q)\ (m 3 + ri)! (m 3 + q)!\ V 2 
V (n — n)! (■ n — q)\ r*i! {q)\ ) 

{n + m 2 - r 2 )! (m 4 + r 2 )! . . 

/ \ 1 1 . c y j 

(n-r 2 )! r 2 ! 

Similarly, one can also find the two party reduced density 
matrix, P 13 { mj }, for photon subtraction by tracing out 
the second and fourth modes in Eq. (14). 

If the first and second modes act as players, we find 
that subtraction is always better than addition for arbi¬ 
trary values of mi and m 2 . This case is similar to the 
case with a single mode being player and cases with sec¬ 
ond and third modes or second and fourth modes being 
players. To show once more that spectators play a fun¬ 
damental role in interchanging the entanglement prop¬ 
erty for photon addition and subtraction, we elaborate 
the analysis in two scenarios - ( 1 ) when a fixed number 
of photons are added (subtracted) in a single spectator 
mode, a positive region emerges, which indicates that the 
quantum correlation in the 13 : 24 bipartition is greater 
for photon addition than that for subtraction, as already 
seen before. An interesting point to note here is that a 
positive region appears for small values of m 2 and almost 
for all values of mi. This is probably due to the fact that 
we add photons in the third mode which belongs to the 
same block as the first mode. (2) When both the specta¬ 
tor modes are active, the positive region can be seen in 
both the axis due to symmetry present in the FM state, 
as depicted in Fig. 10(b). 

Finally, we concentrate on a nontrivial partition, the 
12 : 34 cut (see Fig. 11). From the perspective of en¬ 
tanglement, this partition is unique. In this scenario, 
there are three ways to choose the players. We find that 
with and without participation of spectator modes, en¬ 
tanglement of photon subtraction always higher or equal 


to that of the photon addition which makes this situation 
exclusive from others. 


VI. COMPARISON OF LOGARITHMIC 

NEGATIVITY BETWEEN TWO-MODE AND 
FOUR-MODE STATES 

Upto now, we have considered an FMSV state as in¬ 
put and have compared the behavior of entanglement 
between photon-added and -subtracted states as well as 
entanglement of an output state in different bipartitions 
having different player and spectator modes. In this sec¬ 
tion, our aim is to make comparison between the output 
state obtained from the TMSV state after adding or sub¬ 
tracting photons and the two mode state obtained from 
the FMSV state. To perform such comparison, we dis¬ 
card two modes from the four-mode state and calculate 
the LN of the two mode reduced state, which we then 
compare with the LN of the photon-added (subtracted) 
state that is obtained from the TMSV state as the in¬ 
put [24]. In case of the TM state, the output state, af¬ 
ter adding (subtracting) photons, still remains pure and 
hence LN can be calculated analytically [37]. However, 
for the FM case, the output state is mixed which is ob¬ 
tained by discarding two modes and we adopt the same 
mechanism as we have done to calculate von Neumann 
entropy of reduced density matrices, described in Sec. 
VA3. In particular, we evaluate LN as well as trace for 
large n = N, and then by increasing V, we check whether 
trace goes to unity upto six decimal points. We truncate 
the system when trace has already converged to unity, 
upto six decimal points. 

In the TMSV case, photons can be added to either of 
the modes or to both the modes. On the other hand, 
there are several scenarios for the four-mode states. If 
there is a single player, either one of the mode of the 
output state can act as player or none of the modes of 
the output state is the player. In case of two players, (i) 
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FIG. 12. (Color online) LN betwen the first and second modes 
obtained from the FM output state, (a) First mode is player. 

(b) First and second modes are players with mi + m 2 = 20. 

(c) First mode is player while the second one is spectator with 
m 2 = 10. 


two players can be the two modes of the output state, 
(ii) one mode of the output state can be a player, or (Hi) 
the discarded modes can be the player modes. 

Before considering the FMSV state, let us first con¬ 
sider the TMSV state as input. Note that the nature of 
LN qualitatively matches with the von Neumann entropy 
of the reduced density matrix. As shown in [24], when 
single mode acts as player, the LN for photon addition 
coincide with the subtraction, which is also the case for 
the von Neumann entropy. If both the modes act as play¬ 
ers, photon addition is always beneficial for entanglement 
than the photon subtraction [24]. 

In case of a single player or two players in the FM state, 
if the output state contains the player mode(s), then the 
reduced two-mode state obtained from the photon-added 
state has higher LN than that of the photon-subtracted 
state. Hence, the behavior of LN of the output state from 
TM and FM state are identical. As we have shown, this is 
not the case if we consider the behavior of entanglement 
of pure four-mode output state in bipartitions. Fig. 12 
depicts the behavior of LN of the two-mode reduced state 
from the four-mode output state when the first mode acts 
as player as well as both the modes of the two mode state 
are players. In all these situations, no photons are added 
(subtracted) in the spectator modes. We observe that 
when there is a single player e.g., the first mode of the 
reduced state, entanglement increases (decreases) mono- 
tonically, if photons are added (subtracted). However, 
such monotonicity with respect to the number of pho¬ 
tons added (subtracted) is lost if photons are added (sub- 



FIG. 13. (Color online) (a) The nature of LN of the first 
and second mode against the number of photons added (sub¬ 
tracted) in (from) the third mode which has been traced out. 
(b) LN of the same state in which third as well as fourth 
modes act as players. The equal number of photons are added 
(subtracted) in (from) both the modes. 


tracted) in both the modes with total number of photons 
being fixed as shown in Fig. 12(b). A similar qualitative 
feature in entanglement is seen when the first mode acts 
as player while second mode is a spectator having fixed 
finite number of photons (see Fig. 12(c)). We find that 
the bipartite entanglement reaches its maximum with re¬ 
spect to mi, when equal number of photons are added 
(subtracted) in both the modes i.e. mi = m 2 , in Fig 
12(b) and mi ~ m 2 in Fig. 12(c). 

Lastly, we consider the scenario, when we add and sub¬ 
tract photons in the discarded modes i.e., in the third and 
fourth modes, and we find LN between the first and the 
second modes, which are spectators. LN of the output 
state decreases if one of the discarded modes act as a 
player. For example, by taking third mode as player, 
we plot LN of the first and the second mode with m 3 
in Fig. 13(a). Unlike previous cases, LN of the photon- 
subtracted state is higher than that of the added state 
when m 3 > 9 which can never be observed for the TM 
case. LN of the photon-subtracted state is more pro¬ 
nounced than that of the added one if both the dis¬ 
carded modes act as players. The same number of pho¬ 
tons are added (subtracted) in (from) both the spectator 
modes, i.e. m 3 = 7714 , as shown in Fig. 13(b) in which 
LN(p|f) > LN (plf 1 ). 


VII. NON-CLASSICALITY MEASURE OF THE 
PHOTON-ADDED/-SUBTRACTED FM STATE 


As mentioned in the introduction, the photon addi¬ 
tion and subtraction is one of the ways to create a non- 
Gaussian state. In this section, we quantify the depar¬ 
ture of the photon-added (-subtracted) FMSV state from 
Gaussianity, as a function of added (subtracted) photons 
from the player modes, which was introduced in Sec. III. 

Since the photon-added (-subtracted) FM state is in 
a pure state, the second term of Sng(p): given in Eq. 
(15) vanishes. To calculate Sng(p), we have to find the 
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\ (qm)^z {q2q±)l {qzq^Pz 

where, g* = at + a\, and the expectations are taken over 
the photon-added and -subtracted FM state, given in 
Eqs. (13) and (14) [for details, see the Appendix], 

The Williamson normal form of Eq. (40) can be evalu¬ 
ated by using the prescription given in [38] . We numeri¬ 
cally calculate the Williamson normal form of the matrix 
in Eq. (40) for both photon addition and subtraction 
and calculate the non-Gaussianity, which in this case re¬ 
duces to S(p^^y). In all the cases, photon addition 



m i m-j 


FIG. 14. Behavior of non-Gaussianity measure, 5ng against 
mi. (a) First and second modes are players with mi + m2 = 
10. Spectator modes are ineffective, (b) Spectator modes are 
active with m 3 = 4 and in player modes, mi + m2 = 10. 

leads to a rapid departure of Gaussianity than that of the 
photon subtraction. We also notice that if among four 
modes, photons are added only in two modes, then be¬ 
havior of 6ng obtained in the FM state and the TM state 
are qualitatively similar. It is clear from the behavior 
of the non-Gaussianity measure that photon-subtracted 
state become slowly non-Gaussian as compared to the 
photon-added state and the behavior remains unchanged 
irrespective of the choices of the player and the spectator 
modes (see Fig. 14). The rich picture of the role of differ¬ 
ent modes, captured by entanglement, is not seen by the 
non-Gaussianity measure and hence indicates that there 
is possibly no direct connection between non-Gaussianity 
and entanglement content of the output state obtained 
after photon addition (subtraction) [24]. 

VIII. CONCLUSION 

Photon addition and subtraction constitute useful 
methods to prepare non-Gaussian states. It has already 
been established that non-Gaussian states are useful in 


various quantum mechanical tasks ranging from entan¬ 
glement distillation to quantum error correction. We 
have investigated the entanglement properties of the non- 
Gaussian states generated by adding or subtracting pho¬ 
tons in Gaussian states. In case of two mode states, en¬ 
tanglement of photon-added states are known to be equal 
or higher than that of the photon-subtracted ones. 

We have shown that this is not the case when one in¬ 
creases the number of modes. We found that for four¬ 
mode states, the trend of entanglement distribution in 
different bipartitions of the photon-added (-subtracted) 
states is much richer than that in the two-mode states. 
Specifically, we showed that there exists a scenario, in 
which multimode entanglement content of the photon- 
subtracted state is always higher than that of the corre¬ 
sponding photon-added one. The results remained un¬ 
changed even if one discarded two modes from the four¬ 
mode output state. Moreover, we showed that the pic¬ 
ture that emerges from entanglement of the output state 
does not match with the behavior in the same states of 
distance-based non-Gaussianity measures. Upto now, it 
was known that among addition and subtraction, addi¬ 
tion is more beneficial. But our work shows that photon 
subtraction can also be advantageous if we consider a 
state of a higher number of modes. 
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APPENDIX 

The expectations used in Eq. (40), for the calculation 
of non-Gaussianity, taken over the photon-added and - 
subtracted states are given below. 

(< ll) add = l + 2m 1+ 2 CpMJ 2 (» - ri), (41) 

71,7*1,7*2 
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where E 7l ri r 2 t ^ ie Photon-added states, is the short 
form of E n=0 £"1= Er 2 =o- Further, for the photon- 
subtracted states, 
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